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Abstract 

We calculate the new physics contributions to the rare semileptonic decay B — > X s l + U 
{I = e, u) induced by the charged-Higgs loop diagrams appeared in the top quark two-Higgs 
doublet model (T2HDM). Within the considered parameter space, we found that (a) the effective 
Wilson coefficients C- (m&) (i = 77,91/ and 10 A) in the T2HDM are always standard model 
like; (b) the new physics contributions to and Cgy can be significant in magnitude, but 

they tend to cancel each other; and (c) the T2HDM predictions for Br(B — > X s l + U) agree well 
with the measured value within one standard deviation. 
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I. INTRODUCTION 



Flavor changing neutral current (FCNC) b —>■ s processes are forbidden at the tree 
level in the Standard Model (SM). They proceed at a low rate via penguin or box dia- 
grams. If additional diagrams with non-SM particles contribute to such a decay, their 
amplitudes will interfere with the SM amplitudes and thereby modify the rate as well as 
other properties. This feature makes FCNC processes an ideal place to search for new 
physics. 

In the past decade, the data of B — > X s j decay has served as one of the most impor- 
tant constraints for various new physics models beyond the SM. At present, the world 
average, Br(B X sl ) = (3.55 ± 0.26) x 1(T 4 Q, agrees very well with the standard 
model prediction at next-next-to- leading order (NNLO) 0. The magnitude of the Wilson 
coefficient CV 7 (fib) is therefore strongly constrained by the precision data of B — > X s 
but its sign is still to be determined through the measurement of B — > X s l + l~ decay 
In Ref. |3j, the authors studied B — > XJ + l~ decay and found that the recent experimental 
data of BriB — > X s l + l~) prefer a SM-like CV 7 (/if,). 

In fact, the semileptonic decays B — > X s l + l~ (I = e, /i) have been extensively, in- 
vestigated, for example, in the SM 0,0, the two-Higgs doublet models (2HDM) p] or 
the supersymmetric models 0,13- Our goal in the present work is to calculate the new 
physics contributions to the branching ratio of B — > X s 7 and B — > X s l + l~ decays induced 
by the charged Higgs loop diagrams in the top-quark two-Higgs-doublet model (T2HDM) 
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, and compare the theoretical predictions in the T2HDM with currently available 
data. 

The outline of the paper is as follows. In section II, we give a brief review for the top- 
quark two-Higgs-doublet model and we calculate the new penguin diagrams induced by 
new particles and extract out the new physics parts of the Wilson coefficients or some basic 
functions in the T2HDM. In section III, we present the numerical results of the branching 
ratios of B — > X s l + l~ decay in the SM and the T2HDM, and make phenomenological 
analysis. The conclusions are included in the final section. 



II. THEORETICAL FRAMEWORK 

A. Outline of the top quark two-Higgs-doublet model 

The specific model considered here is the top quark two-Higgs-doublet model (T2HDM) 
proposed in Ref. and studied in Refs. 1(], U|, which is also a special case of the 2HDM 



of type III In this model, the large mass of the top quark arises naturally in the 

extension of the SM since the top quark is the only fermion receiving its mass from the 
vacuum expectation value (VEV) of the second Higgs doublet. All the other fermions 
receive their masses from the VEV of the first Higgs doublet. 

Let us now briefly recapitulate some important features of the model of Ref. 0]. Con- 
sider the Yukawa Lagrangian of the form: 

C Y = -L L faEl R - Q L <piFd R - Q L faGlV*u R - Q L 4> 2 Gl {2) u R + H.c. (1) 

where Ql and Ll are 3- vector of the left-handed quark and lepton doublets, respectively; 
4>i {i = 1, 2) are the two Higgs doublets with (pi = iT 2 (j)*; and E, F,G are the 3x3 matrices 
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in the generation space and give masses respectively to the charged leptons, the down and 
up type quarks; = diag(l, 1, 0); = diag(0, 0, 1) are the two orthogonal projection 
operators onto the first two and the third families respectively The top quark is assigned 
a special status by coupling it to one Higgs doublet that gets a large VEV, whereas all 
the other quarks are coupled only to the other Higgs doublet whose VEV is much smaller. 
Consequently, if one sets the VEVs of fa and 2 to be v x j\[2 and v 2 e w /V^ B, the ratio 
of two Higgs VEVs, tan/3 =Vifv\, is required to be relatively large. 

The Yukawa couplings involving the charged-Higgs bosons are of the form 

= -jM- — {-u L VM D d R [G + - t&n(3H + ] + u R MuVd L [G + - tan/3# + ] 

+u R E f Vd L [t&ni3 + cot (3}H + + h.c.}. (2) 

where G ± and H ± denote the would-be Goldstone bosons and the physical charged Higgs 
bosons, respectively. Here Mu and Mp are the diagonal up- and down-type mass matrices, 
V is the usual CKM matrix and S = MtjU R l^U R . U R is the unitary matrix which 
diagonalizes the right-handed up-type quarks and has the following form: 



U R = sin0 cos0 0x0 V 1 - \zM 2 -ectC ■ (3) 





where e ct = m c /m t , £ = \^\e tS is a complex number of order unity, and the phase 5 in £ is 
a new CP violating phase. Inserting Eq. (J5J) into the definition of S yields 





m c el\i\ 2 m^V 1 - I^CI 2 I • (4) 
m c ^l - |e ci e| 2 m t {l - |e ct £| 2 ) 

In the following sections, we will calculate the charged Higgs contributions to the rare 
decay B — > X s l + l~ in the top quark two-Higgs-doublet model. 



B. Effective Hamiltonian for B X s l+l~ in the SM 

In the framework of the SM, the effective hamiltonian inducing the transition b — > sl + l~ 
at the scale \i can be written as follows: 

AC 10 

n = (5) 
* i=i 
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where Gf is the coupling constant, and V t *Vtb is the Cabibbo-Kobayashi-Maskawa (CKM) 
factor The operators can be chosen as Ref. jij: 

Qi = (^7^ a CL)(cL7 At T a 6 L ) , Q 2 = (s Lltl c L )(c L Yb L ) , 

q i 

Qi 1 = -^{s^b^F^ , Q 8g = - mb (s L a^T a b R )G% , 

9 S 9s 

2 2 

Qw = - 2 {sliM 5>~^) ' Git* = - 2 ^Li,b L ) Ytf-f-rsf) > (6) 

9 S t 9 S e 

where Qi i2 are the current-current operators, Q3-6 the QCD penguin operators, Q? i8 
"magnetic penguin" operators, and Q9,io semileptonic electroweak penguin operators. 
T a (a = 1,...,8) stands for SU(3) C generators, L,R = (1 q= ^y 5 ) / 2 by definition. The 
sum over q runs over the quark fields that are active at the scale \i = 0(m,b), i.e., 
q G {u, d, s, c, b}. We work in the approximation where the combination (V* s V u b) of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements is neglected. We do not separate 
top-quark and charm-quark contributions and will give the results in the summed form. 

To calculate the semileptonic B meson decays at next-to- leading order in a s , we should 
determinate the Wilson coefficient Ci(Mw) through matching of the full theory onto 
the five-quark low energy effective theory where the W gauge boson, top quark and 
the new particles of T2HDM heavier than M\y are integrated out, and run the Wilson 
coefficients down to the low energy scale \i ~ Oiviib) by using the QCD renormalization 
group equations. The corresponding Wilson coefficients in SM can be found, for example, 
in Refs.jllQ. 



C. New physics contributions 

In the framework of the SM, the semileptonic B — > X s l + l~ (/ = e~,/i~) decays pro- 
ceed through loop diagrams and are of forth order in the weak coupling. The dominant 
contributions to this decay come from the W box and Z penguin diagrams. The corre- 
spondin g on e-loop diagrams in the SM were evaluated long time ago and can be found in 
Refs. pL lKjl]. The calculations at the next-next-to-leading order (NNLO)are also available 
now. 

In the T2HDM considered here, besides the SM diagrams with a W-gauge boson and 
an up quark in the loop, the B — > X s l + l~ decays can also proceed via the new diagrams 
involving the charged-Higgs boson exchanges, as illustrated by Fig^ In order to deter- 
mine the new physics contributions to the relevant Wilson coefficients C 77 , C 8g , C$v, and 
Ciqa at the Mw scale, we need to calculate the corresponding Feynman diagrams. 

The new physics parts of the Wilson coefficients CV 7 and C 8g have been calculated in 



Refs. jlOL and confirmed by our independent calculation. In the naive dimensional 
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•7,Z 



>W 



i i i n i i v x i 

(d) (e) (f) 

FIG. 1: The typical Feynman diagrams for the decay B — > X s / + Z~ in the T2HDM. The internal 
solid and dashed lines denote the propagators of upper quarks (u,c,t) and charged Higgs boson, 
respectively. 



regularization (NDR) scheme, they are of the form 

1 



C^(M W ) 



i=c,t 



tan 2 (3 + — — (E 7 V%(tan 2 /3 + 1) 



1 



5(1/,) + -Afo) 
o 



miV is 



■1 + 



1 



-(EVUcot^+l) 



i=c,t 



- tan 2 /3 + 



E( Vi ) + -F( Vi ) 
o 



-1 + 



(S V%(tan 2 /3 + l) 



■(EV), b (cot 2 /? + l) 



(6 



with the Inami-Lim functions 

= 

^(i/) = 
5(y) = 



7y-5y 2 -8y 3 2y 2 - 3y 3 

12(1 -y)3 + 
-3y + 5y 2 2y - 3y 2 



12(1 -y) 2 6(1 -y) 3 
-3y + y 2 



Hv], 



y 

4(1 -yy 2(1 -y)3 



My], 



+ 



3^ 



% + 5y 2 - 
4(1 -y) 3 ' 2(1 -y) 4 



ln[y], 



(7) 



(9) 



where K is = -V lb V*J(V tb V t * s ), yi = (mi/m H ) 2 . 

As for the Wilson coefficients Cqv, and CW at the Mw scale, we found the new physics 
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parts after calculating analytically the Feynman diagrams as shown in Fig. ^ 

1 



w 



sin 6 W 
1 



[C^-Bf]-[Df + 4C l 
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sin^fl 
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rnimj, tan 2 (5 



E 
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„ mf 



K 



i=c,t 



8M 2 , 



B + (x H +,x t ), 
C oi(Vi) 



(10) 
(11) 

(12) 



4m 2 , . o n n' / \ 

VwC u {yi) 



- tan 2 p + 



-1 + 



2 

mf 



(1 - 3 sin 2 ^)C 01 (yi) - C n (y t ) 



- tan 2 /3 + — — (E T V%(tan 2 /? + 1) 



(13) 



,2^(2/0 



i=c,t 



- tan 2 /3 + — L-(£ T y%(tan 2 (3 + 1) 
miV is 



-1 + 



-(£VWcot 2 /? + l) 



(14) 



B+{x,z) 

H(y) 
c'oM 

c n (y) 



x — z 



\n[z] ln[x] 



z — 1 a; 



38y - 79y 2 + 47y 3 4y - 6y 2 + 3y 4 



72(1 -y)3 12(1 -y)* 



+ 



1-y (1-y) 2 
3y - y 2 y 



ln[y], 



+ 



4(1 -y) 2 2(1 -y)3 



(15) 



where yj = m 2 /m 2 H , xh+ = mfj/M^, and x t = m 2 /M^. V is the CKM matrix, and the 
matrix E has been given in Eq. (J3J . The contributions from Fig^ and the Fig^ when 
the internal W and charged-Higgs lines exchange their position are strongly suppressed 
by a factor of (mi/mn) 2 {m\ = m e ,m M ) or m s /mb, and therefore have been neglected. 
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D. The differential decay rate 



Within the Standard Model, the differential decay rate for the decay B — > XJ + l in 
the NNLO approximation can be written as jH 0] 



R(s) = dS 



4r(6 - sin- 



r(6 — > cez/ 
2 



+4 1 + 



Oy 7 



47T 2 
2 



12Re 



(1-g) 2 

f(z)K(z) 



(l + 2s) 



^ioa(^) 



7 7 I'-w I s 



where 



that are related to the evolved coefficients Ck(fib) as follows: 

20 



with 
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C e 9 f v f (s) 



4tt 14 
-——C 7 {fi b ) - -C 3 (n b ) - -C 4 (/i b ) - —C 5 (ii b ) 
a s [ii b ) 3 9 3 



80 
~9~ 



4C 9 (/x 6 ) 



7T 



-CM + C 2 (ji b ) + 6C?(/i 6 ) + 60C 5 Q (/i fe ) 



+h — §,S 
W 

/ 7 9 32 

+7i(M) (--C 3 (/i6) - -38C 5 (/i 6 ) - yC 6 (/i 6 ) 



+/i(0,s) 



1 2 

2^3(^6) - ^C 4 (fi b ) - 8C 5 (/i fc ) 
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4 64 64 

+ 3^3(^6) + -gC 5 (fi b ) + — C 6 (/x 6 ), 



Cioa{s) 



\fjeff 



4Cio(^fe) 



7T 



m. 



+ uj{s) 



h(0,s)-h[ -§,s 



7?>r 
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4 

-Ci(/i 6 ) + C 2 (/i6) 



--(2 + x)v/|l-x| 



In 



itc. 



for x = 4z/s < 1, 



2 arctan(l/\/x — 1), for x = 4z/s > 1, 



h{0,s) 



8 4 n - ■ ^ 
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4 r . 2, 2 2 5 + 4s 
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5 + 9s - 6s 2 
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C 6 (/x 6 ), (18) 
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(20) 

(21) 
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2a s (/i) 
3tt 



7T 



- 2Az 4 In 2 
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(1 



(25) 
(26) 



here s = (pi+ + pi-) 2 /ml = mjjm 2 ., z = m c /m b , f(z) and k(z) are the phase-factor and 
single gluonQCD correction to the b — > ceu decay, respectively. 



In Refs. QJ 



the Wilson coefficients have been expanded perturbatively as follows 
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C; 



d 0) + 
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9* 



(An) 



CP + 0(g 



6^ 



(27) 



For the standard model parts of the Wilson coefficients C^^Cf'' and Cf 1 the explicit 
expressions as given in Refs. 0, 0| will be used in our numerical calculation. For 
the new physics part, only C\ (Mw) are known at present, as given explicitly in 
Eqs. (l7l8ll(Jlllj) . and will be included in numerical calculations. 



t(2) 



III. NUMERICAL RESULT 

In this section, we first give the input parameters needed in numerical calculations, 
and then present the numerical results and make some theoretical analysis. 



A. input parameters 



In numerical calculations we will use the following input parameters (all masses are in 
GeV) H3: 



M w = 80.425, Gp = 1.16639 x l(T 5 GeV 



-2 



a,. 



1/128, 



A® =0.22r, 



m c = 1.4, m 6 = 4.8±0.2, m t = 173.8 ± 5, 
A = 0.853, A = 0.2200, p = 0.20 ± 0.09, rj = 0.33 ± 0.05, 
sin 2 ^ = 0.23124, Br(B -> X c ev) = 0.1061, 



(28) 



where the parameter A, A, p and r] are Wolfenstein parameters of the CKM mixing matrix. 
For the strong coupling constant a s (p) we use the two-loop expression, 



a s (p) 



An 



with 



0o 



/% ln(/i7A^ 
33-2/ 



lnhuy/A. 



P 2 o ln(/x 2 /A 



—) 

US' 



(29) 



Pi = 72- 10 / - 8//3 



(30) 



where the / is the number of quark flavors and the term MS denotes the modified sub- 
traction scheme. 
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B. B — > Xg'j decay 

There are four free parameters m#, tan/5, |£| and a new CP-violating phase 5 in 
the T2HDM. We fix |£| = 1 throughout the paper and consider other three as variable 
parameters to be constrained by precise measurements, such as the date of Br(B — > X s ^). 

In Ref. [HI, the branching ratio Br(B — > X s j) have been calculated in both the SM 
and the T2HDM. Using the formulas as given in Appendix El and taking the range of 

2.77 x 10" 4 < Br(B -> X sl ) < 4.33 x KT 4 (31) 

as the experimentally allowed region at 3a level one can read off the lower limit on 
the mass of charged-Higgs boson m# directly from Fig. EJ 

m H > 300GeV, (32) 

for fixed tan j3 = 30 and 5 = 0°. It is easy to see from Fig. |21 that (a) a light charged Higgs 
boson with a mass less than 200 GeV is excluded by the data of B — > X s j decay at 3a 
level; and (b) a charged-Higgs boson with a mass heavier than 300 GeV is still allowed 
by the same data. 

6.0 i i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i 

5.5 1 

. „ *V\\ \ 10 : 

5.0 -:■ , ■. 

:\\ \ •.2o'--.. : 




2 f .... i .... i .... i .... i .... i .... i .... i .... ' 

200 300 400 500 600 700 800 900 1000 

m H (GeV) 



FIG. 2: The m H dependence of Br(B X s -f) in the T2HDM for 5 = 0°, and for tan /3 = 10, 
20, 30, 40 and 50 respectively. The band between two horizontal dash dot lines shows data as 
specified in Ea. H31|) . The solid horizontal line shows the central value of the SM prediction. 

As shown in the contour plot Fig. |HJ the region between the short-dashed and solid 
curves is still allowed by the data of B — > Xg'j as given in Eq. (J3~T|) for fixed value of 
5 = 0°. On the other hand, by assuming tan/5 = 30 and ran = 400 GeV, one finds a 
strong constraint on the phase 5: 5 < 44°. 
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FIG. 3: Contour plot in tan (3 — mu plane obtained by considering the data in Eq. I)31j) for fixed 
5 = 0. The region between the short-dashed and solid curves is still allowed by the data as given 
in Eq. EU 



C. B — > X s l+r decay 
The branching ratio of B — > X s l + l~ (I = e, /i) has been recently measured by BaBar 



and Belle Collaborations[20, 2lj. In the low-g 2 region 1 , the average of BaBar and Belle's 
measurements is 

Br(B -> X s / + r) = (1.60 ± 0.51) x 10~ 6 . (33) 
Theoretically, the integrated branching ratio can be written as 

Br u = Br(B -> XJu) ! " R(s), (34) 

where s = q 2 /ml with s a = I /ml and Sf, = 6 /ml, and the differential decay rate R(s) 
has been defined in Eq. (jib)) . The SM prediction after integrating over the low-g 2 region 
reads 

Br{B -> X s /+r) = (1.58 ± 0.08| mi ± 0.07| Mb ± 0.04| Cim ± 0.06| mfc + 0.18| M J x lO" 6 
= (1.58 ± 0.13 + 0.18| MW .) x HT 6 . (35) 

where the errors correspond to the uncertainty of input parameters of nit, A P, V and 
nib as shown in Eq. (}2"Hj) . and for m^jl < fib < 2m&. The last error refers to the choice 
of 

pw = 120 GeV, instead of pw = Mw- Since we here focus on the new physics 



1 The low-q 2 region is the region with lGeV 2 < = q 2 < 6GeV 2 . 
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TABLE I: The effective Wilson coefficients and the interference term (12Re[C^ (Cgy )*]) for 
fixed s = q 2 /m 2 = 0.2, the branching ratio integrated over the low-q 2 region in units in 
the SM and the T2HDM for m H = 300, tan/3 = 10,30,50 and <5 = 0° (a), 30° (b) and 60° (c). 
Only the central values are shown here. 







r*ff 

a V 


10A 


Int. Term 


Bru 


SM 


-0.344 


4.302 + £0.064 


-3.547 


-17.73 


1.579 


T2HDM 

tan (5 = 10 


(a) -0.422 + i0.001 

(b) -0.424 + £0.006 

(c) -0.428 + £0.010 


(a) 4.205 + £0.063 

(b) 4.218 + £0.014 

(c) 4.255 - £0.021 


(a) -3.552 

(b) -3.552 + £0.001 

(c) -3.553 + £0.001 


-21.30 
-21.45 
-21.84 


1.576 
1.581 
1.595 


tan/? = 30 


(a) -0.376 + £0.001 

(b) -0.389 + £0.050 

(c) -0.425 + £0.086 


(a) 3.430 + £0.051 

(b) 3.554 - £0.385 

(c) 3.879 - £0.700 


(a) -3.546 

(b) -3.546 + £0.001 

(c) -3.547 + £0.002 


-15.47 
-16.84 
-20.50 


1.342 
1.388 
1.581 


tan/? = 50 


(a) -0.283 + £0.002 

(b) -0.321 + £0.140 

(c) -0.420 + £0.239 


(a) 1.882 + £0.026 

(b) 2.226 - il.183 

(c) 3.126 - £2.056 


(a) -3.544 

(b) -3.544 + £0.002 

(c) -3.546 + £0.003 


-6.40 

-10.56 

-21.65 


1.033 
1.167 
1.526 



contributions to the branching ratios of B — > X s l + l~ decay, we will take nw — m 
the following without further specification. 

Now we consider the new physics contributions. When the new physics parts of the 
Wilson coefficients cf ){Mw) for i = 77, 8g, 9V and 10 A are taken into account, the values 
of the effective Wilson coefficients appeared in Eq. (fTrjj) and the theoretical predictions of 
the branching ratio will be changed accordingly, as listed in Table |l]for tan/3 = 10, 30, 50, 
m H = 300 GeV and 5 = 0°, 30° and 60° . 

In Figs. 0] and O in order to show more details of the m# and tan f3 dependence, we 
draw the real part of the effective Wilson coefficients {nib) and Cgy (mb) for fixed 
s = q 2 /ml = 0.2 and 5 = 0°. Within the considered parameter space of the T2HDM, it 
is easy to see from the numerical results in Table H] and Figs. 0] and El that 

(i) The effective Wilson coefficient C^(m&) is always SM-like. This feature can be 

seen explicitly in Fig. HI where the m#-dependence of the real part of Cj (m&) is 
shown for 5 = 0°, tan/5 = 10, 30, 50 and 200GeV <m H < lOOOGeV. The imaginary 
part of (nib) is generally small. 

(ii) The effective Wilson coefficient Cgy (nib) is also SM-like. The imaginary part 
of Cgy(mb) is also generally small. 

(iii) The new physics contribution to C^ A is very small in size, less than 1% of its 
standard model counterpart, and therefore can be neglected safely. 

(iv) The new physics contributions to and Cgy can be significant in magnitude 
respectively for large tan f3, large S and lighter charged-Higgs boson, as can be seen 
from the numerical results in Table Eland illustrated explicitly by FigsIUandEl But 
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FIG. 4: The m# dependence of the real part of the effective Wilson coefficient (ra^) in the 
SM (solid line) and T2HDM for 5 = 0°, and tan/3 = 10 (dots curve), 30 (dot-dashed curve) and 
50 (dashed curve), respectively. 




200 300 400 500 600 700 800 900 1000 

M H (GeV) 



FIG. 5: The mu dependence of the real part of the effective Wilson coefficient C^y {m^) in the 
SM (solid line) and T2HDM for s = 0.2, 5 = 0°, and tan/3 = 10 (dots curve), 30 (dot-dashed 
curve) and 50 (dashed curve), respectively. 

they tend to cancel each other and finally lead to a small change to the prediction 
for the branching ratio under study. 

It is worth noting that both the real and imaginary parts of effective Wilson coefficients 
are taken into account in our calculation of the branching ratio. 

As shown in Eq. f!16|) . the differential decay rate depends on the summation of three 
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FIG. 6: The run dependence of the branching ratio of B — > X s l + l~ in the SM and T2HDM for 
5 = 0°, and tan/3 = 30. The contributions from the term-1, term-2, interference term and their 
summation are shown by the dot-dashed, dashed, short-dashed and solid curve, respectively. The 
horizontal band between two dots line shows the data: Br(B — » X s l + l~) = (1.60 ±0.51) x 10 -6 , 
while the solid line refers to the central value of SM prediction: Br{B — > X s l + l~) = 1.58 x 10 -6 . 

terms: 



Term — 1 : 

Term — 2 : 
Term — 3 : 




where the third term is the interference term, which has opposite sign compared to first 
two terms. From Fig. El one can see easily that 

(i) After the inclusion of new physics contributions in T2HDM, the signs of three 
terms remain unchanged. 

(ii) The new physics contributions to these three terms are indeed tend to cancel 
each other and result in a summation (solid curve in Fig. |UJ) which becomes closer to 
the SM prediction (solid line in Fig. IHJ) when mn becoming larger. The theoretical 
predictions for the branching ratio in the SM and T2HDM agree well for the whole 
range of mn considered here. They are also in good agreement with the data within 
one standard deviation. 

Analogous to Fig. El the Figs. [7| and |S] show the tan/3 and S— dependence of the branch- 
ing ratio Br(B — > X s l + l~), respectively. Here, the cancelation of new physics contribu- 
tions to different terms occurs and leaves the summation, the theoretical prediction in 
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FIG. 7: The same as Fig. H3 but shows the tan j3 dependence of the branching ratio of B 
X s l + l- in the SM and T2HDM for S = 0° and m H = 300 GeV. 
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FIG. 8: The same as Fig. H3 but shows the S dependence of the branching ratio of B — > X s l + l 
in the SM and T2HDM for tan (5 = 30 and m H = 300 GeV. 



the T2HDM, in good agreement with the SM prediction as well as the measured value 
within one standard deviation. From Fig. [7J one can also see that a tan (3 smaller than 
40 is preferred by current data. 
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IV. SUMMARY 



In this paper, we calculate the new physics contributions to the branching ratio of 
B — > X s 7 and B — > X s l + l~ decays induced by the charged Higgs loop diagrams in 
the top-quark two-Higgs-doublet model, and compare the theoretical predictions in the 
T2HDM with currently available data. 

In Sec. |nj we firstly present a brief review about the basic structure of the top-quark 
two-Higgs-doublet model, and then evaluate analytically the new Feynman diagrams in- 
duced by the charged Higgs H ± exchanges and extract the new physics parts of the 
Wilson coefficients C^(/j,w), Csgif^w), Cgv(fJ'w) an d Cwa^w) which govern the new 
physics contributions to B — > X s 7 and B — > X s l + l~ decays considered in this paper. For 
the SM part, we use the known analytical formulae at NNLO level as given for example 
m Refs. |5j,llJ. The new physics contributions are included through the modifications to 
the corresponding Wilson coefficients at matching scale nw ~ Mw- 

From the numerical results and the figures as shown in Sec. Hill we found that 

(i) For the T2HDM studied here, a light charged Higgs boson with a mass less than 
200 GeV is excluded by the data of B — > X s ^ decay at 3cr level. But a charged- 
Higgs boson with a mass heavier than 300 GeV is still allowed by the data of both 
B — > Xs'y and B — > X s l + l~ decay. The data of B — > X s ^ also prefer a small 5, a 
new CP violating phase appeared in the Yukawa couplings of the T2HDM. 

(ii) After the inclusion of new physics contributions, the effective Wilson coefficients 
C\ (nib) (i = 77,9V and 10A), which govern the branching ratio of B — > X s l + l~ 
decay, are always SM-like within the considered parameter space of T2HDM. The 
sign of the interference term in Eq. (jl6j) remains unchanged. 

(iii) The new physics contributions to and Cg{f can be significant in magnitude 
respectively for large tan (3, large S and lighter charged-Higgs boson, but they tend 
to cancel each other and finally result in only a small change to the prediction for 
the branching ratio of B — > X s l + l~ decay. This feature can be seen clearly through 
the numerical results in Table and the curves shown in last three figures. 

(iv) Within the considered parameter space of the T2HDM, the T2HDM predictions 
for BriB — > X s l + l~) agree well with the SM as well as the measured value within 
one standard deviation. 
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APPENDIX A: Br(B X s -y) IN THE SM AND T2HDM 



The branching ratio of B — > X s ^y at the next-to-leading order (NLO) in the SM and 
the leading order (LO) in the T2HDM can be written as |igL liT 



B(B - X sl ) = E S L 



v; s v tb 



cb 



6a P 



nf(z)K(z) 



\D\ 2 + A + A] 



(Al) 



where Esl = 10.61% is the measured semileptonic branching ratio of B meson, a em =l/128 
is the fine-structure constant, z = mP° le / 'm p h ole = 0.29 ± 0.02 is the ratio of the quark pole 
mass. The function f(z) and k(z) have been given in Eqs. (|25|) and ()26jl. 

The term D at low energy scale \x = 0(rrib) in Eq. ()A1|) corresponds to the subprocess 
b — ► S7 



(A2) 



Here C™(n) denotes the SM part of the Wilson coefficient CV 7 (/i) at NLO level, and the 
explicit expression of C^(ji) at both LO and NLO level can be found easily in Ref. (l5l |. 

The new physics part of the Wilson coefficient CV 7 and Cs g at the matching scale Mw 
are currently known at LO level and have been given in Eqs. (JJJ) and (JHJ). At the low 
energy scale \i = 0(rrib), the leading order Wilson coefficients C^{p) and Cg 5 P (/z) can be 
written as 

C-Gu) = #C-(M ff ) + ~ - r/I) C 8 N ;(M W ), (A3) 

where rj = a s (Mw)/cx s (fj), and the Wilson coefficient C^(Mw) has been given in Eq. (jSJ). 
The function V"(/i) in Eq. ()A1|) is defined as 2^1 



r(//) = ^£ C o (/1 ) 



47T 



i=l 



r i+2 7 j7 ln^ 2 



(A5) 



where the functions r, (i = 1, . . . , 8) are the virtual correction functions (see Appendix D 
of Ref. |22|), 7 4 °7 are the elements of the anomalous dimension matrix which govern the 
evolution of the Wilson coefficients from the matching scale Mw to lower scale \x. The 
values of 7° 7 can be found in Ref. [22l ]. 

In Eq. ()Al|) . the term A = A(fi) is the the correction coming from the bremsstrahlung 
process b — > sjg 23| 



a s (n) 



7T 



'■J 



Re{^°(/i) [C°(/i)]* /y} . 



(A6) 



The coefficients have been defined and computed in Refs. [2 
explicit expressions of those relevant /y as given in Appendix E of Ref. 2^ . 
Finally, the term A in Eq. (|A1J) denotes the non-perturbative corrections 12 



3. We here use the 



5 NP 

7 

ml 



(A7) 
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with 

S? = J-l x *> C = -y, (A8) 

where A 2 = (m|» - m|)/4 = 0.12 GeV 2 and Ai = 0.5 GeV 2 . 

In the expressions of V(/x), A(fi) and A, the superscript "0" means that the correspond- 
ing Wilson coefficients at LO level will be used. The numerical results show that the new 
physics contributions to "small quantities" A(fi) and A are very small in magnitude and 
can be neglected safely. 
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